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We describe an experiment that has set new limits on the time reversal invariance violating D 
coefficient in neutron beta-decay. The emiT experiment measured the angular correlation (J) ■ (p e x 
p p ) using an octagonal symmetry that optimizes electron-proton coincidence rates. The result is 
D — [— 0.6±1.2(stat)±0.5(syst)] x 10 -3 . This improves constraints on the phase of gA/gv and limits 
contributions to T violation due to leptoquarks. This paper presents details of the experiment, data 
analysis, and the investigation of systematic effects. 

PACS number(s): 11.30.Er, 12.15.Ji, 13.30. Ce, 23.20 En 



I. INTRODUCTION 

CP violation has been observed so far only in the de- 
cays of neutral kaons JjJ]. Recently evidence for the im- 
plied T violation in the neutral kaon system has been re- 
ported [Q. These effects could be due to the Kobayashi- 
Maskawa phase in the Standard Model ||. However, 
these observations could also be due to new physics, and 
it is well-established that new sources of CP violation are 
required by the observed baryon asymmetry of the uni- 
verse Many extensions of the Standard Model con- 
tain new sources of CP violation and can be probed in 
observables for which the contribution of the Kobayashi- 
Maskawa phase in the Standard Model is small. The 
present experiment searches for CP violation in one such 
observable, a T-odd correlation in the decay of free neu- 
trons. 

The differential decay rate for a free neutron can be 
written |Hl 
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conservation of momentum allows the transformation of 
the triple-correlation term into 
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where p p is the momentum of the recoil proton. 

The D coefficient is sensitive only to T-odd interactions 
with vector and axial vector currents. In a theory with 
such currents, the coefficients of the correlations depend 
on the magnitude and phase of A = |A|e _t ^, where |A| = 
\9a/9v\ is the magnitude the ratio of the axial vector to 
vector form factors of the nucleon. In this notation, the 
coefficients are given by 
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The most accurate determinations of |A| (current world 
average |A| = 1.2670 ± 0.0035) come from measurements 
of A B. The coefficients a, A, and B, respectively, are 



where p e , E e and p v , E v are the momentum and energy 
of the outgoing electron and neutrino, respectively, S(E e ) 
is a phase space factor, and (J) is the neutron spin. The 
triple-correlation D(3) ■ (p e x p„) is odd under motion 
reversal, and can be used to measure time reversal in- 
variance violation when final state interactions are taken 
into account. Note that in the rest frame of the neutron, 
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measured to be -0.102 ± 0.005, -0.1162 ± 0.0013, and 
0.983 ± 0.004 ||. Several previous experiments found the 
value of -D, and thus sine/), to be consistent with zero at 
a level of precision well below 1%. The three most recent 
such measurements found D = (— 1.1±1.7) x 10~ 3 and 
D = (2.2 ± 3.0) x 10~ 3 ||, and D = (-2.7 ± 5.0) x 10" 3 
@, constraining <j) to 180.07° ± 0.18° §. 

Final state interactions give rise to phase shifts of the 
outgoing electron and proton Coulomb waves that are 
time reversal invariant but motion reversal non- invariant. 
Thus D has terms that arise from phase shifts due to pure 
Coulomb and weak magnetism scattering. The Coulomb 
term vanishes in lowest order in V- A theory |5 1 , but scalar 
and tensor interactions could contribute. The Fierz in- 
terference coefficient measurements (ltjjllj can be used 
in limiting this possible contribution to 
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Interference between Coloumb scattering amplitudes and 
the weak magnetism amplitudes produces a final state 
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The D coefficient has also been measured for 19 Ne de- 
cay, with the most precise experiment finding D^e = 
(4 ± 8) x 10~ 4 |0|. The predicted final state effects 
for 1 Ne are approximately an order of magnitude larger 
than those for the neutron and may be measured in the 
next generation of 19 Ne experiments. For 8 Li, a triple- 
correlation of nuclear spin, electron spin and electron 
momentum has been measured, with the most precise 
measurement at R = (0.9 ± 2.2) x 10~ 3 |]|. Unlike 
D, a nonzero R requires the presence of scalar or tensor 
couplings and thus is a tool to search for such couplings. 
The electric dipole moments (EDMs) of the electron [[l5| , 
neutron [0, and 199 Hg atom @ are arguably the most 
precisely-measured T-violating parameters and bear on 
many of the same theories as D. Table [| summarizes 
the current constraints on D from analyses of data on 
other T-odd observables for the Standard Model and ex- 
tensions ]is|] . For lines 2-5 these limits are derived from 
the measured neutron or 199 Hg EDM. 
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1. Kobayashi-Maskawa Phase 

2. Theta-QCD 

3. Supersymmetry 

4. Left-Right Symmetry 

5. Exotic Fermion 

6. Leptoquark 
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< 10" 14 
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TABLE I. Constraints on D from analyses of other T-odd 
observables for the Standard Model and extensions. 



In the nearly two orders of magnitude between the 
present limit on D and the final state effects lies the op- 
portunity to directly observe or limit new physics. More- 
over, accurate calculations of magnitude and energy de- 
pendence of the final state effects can be made to extend 
the range of exploration still further. 



II. OVERVIEW OF THE EMIT DETECTOR 

In the emiT apparatus, a beam of cold neutrons is 
polarized and collimated before it passes through a de- 
tection chamber with electron and proton detectors (four 
each). A schematic of the experiment is shown in Figure 
|l|. The most significant improvements over previous ex- 
periments are the achievement of near-unity polarization 
(> 93% compared to 70% in 0) and the construction of 
a detector with greater acceptance and greater sensitivity 
to the D coefficient. The octagonal arrangement of the 
eight detector segments gives them nearly full coverage 
of the 27T of azimuthal angle around the beam, nearly 
twice the angular acceptance in previous experiments, 
and the detector segments are longer than in previous 
experiments. The placement of the two types of detec- 
tors at relative angles of 135° is also an improvement over 
previous experiments, in which the coincidences were de- 
tected at 90°. While the cross product is greatest at 90°, 
the preference for larger electron-proton angles in the de- 
cay makes placement of the detectors at 135° the best 
choice to achieve greater symmetry, greater acceptance, 
and greater sensitivity to D (see Figure 0). Combined 
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FIG. 2. Although the cross-product (dashed line) is max- 
imized at electron- proton detection angles of 90°, the overall 
sensitivity to D (solid line) is enhanced at larger angles due to 
the phase space for the decay. Placing the detectors at 135° 
allows for an octagonal geometry that combines greater sym- 
metry, acceptance, and sensitivity when compared to place- 
ment of the detectors at 90°. The solid curve in this figure 
is the sensitivity for a zero-radius beam, which would exhibit 
a factor of 7 enhancement for 135° as compared to 90°. For 
our nearly 3 cm-radius beam, the enhancement factor is close 
to 3. 
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FIG. 1. The emiT experimental apparatus beamline layout. The neutrons traveled through 8.8 meters of guides and vacuum 
components before reaching the beam stop. 



with the higher neutron polarization from the supermir- 
ror polarizer our geometry provides for an overall sen- 
sitivity to D that is a factor of w 7 greater than previ- 
ous measurements, assuming the same cold neutron beam 
flux. 

The first run of the experiment was conducted at the 
NIST Center for Neutron Research (NCNR) in Gaithers- 
burg, MD. The experimental apparatus is outlined below, 
while more detailed descriptions can be found in |]l9| , po| . 

A. Polarized Neutron Beam 

The NCNR operates a 20-MW, heavy-water- 
moderated research reactor. Neutrons from the reactor 
pass through a liquid hydrogen moderator to make cold 
neutrons with an approximately Maxwellian velocity dis- 
tribution at a temperature of about 40 K. The average 
neutron velocity is about 800 m/s. The neutrons are 
transported 68 meters to the apparatus via a 58 Ni-lined 
neutron guide. Neutrons are totally internally reflected 
if they enter with an angle of incidence less than 2 mrad 
for each A of de Broglie wavelength. The capture flux 
of the neutrons was measured using a gold foil activation 
technique to be p n un = 1.4 x 10 9 n/cm 2 /s (where vo = 
2200 m/s) at the end of the neutron guide. (The capture 
flux quantifies the neutron density in the detector for the 
polychromatic beam.) The beam passes through a cryo- 
genic beam filter of 10-15 cm of single crystal bismuth 
which filters out residual fast neutrons and gamma-rays. 

The neutrons are polarized with a double-sided bender- 
type supermirror polarizer obtained from the Institut 
Laue-Langevin in Grenoble, France |2lJ]. The supermir- 
ror consists of 40 Pyrex |^2j plates coated on both sides 
with cobalt, titanium, and gadolinium layers that maxi- 
mize the reflection of neutrons with the desired spin state 
while absorbing nearly all neutrons of the opposite spin 
state. The supermirror was measured to polarize a 4.5 



cm by 5.5 cm beam with 24% transmission relative to the 
incident unpolarized flux. The neutron polarization was 
determined to be > 93% (95% CL). 

The neutrons travel the one meter from the polarizer 
to the spin-flipper inside a Be-coated glass flight tube in 
which a small helium overpressure is maintained to mini- 
mize beam attenuation via air scattering. The neutrons, 
which have spins that are transverse to their motion, then 
pass through two layers of aluminum wires which com- 
prise the current-sheet spin flipper. When the current in 
the second layer is antiparallel to that in the first there 
is no net magnetic field and the neutron polarization is 
unaffected. When the currents are parallel, the neutron 
spin does not adiabatically follow the rapid change in 
field orientation and thus the sense of (J) • B is reversed. 
Downstream of the spin flipper, weak magnetic fields adi- 
abatically rotate the spin to longitudinal, i. e. parallel or 
antiparallel to the neutron momentum. The longitudi- 
nal guide fields are 2.5 mT upstream and 0.5 mT within 
the detector. Figure |^ shows the spin transport system. 
The polarization direction is reversed every 5 seconds. In 
the detection region, the longitudinal field is produced by 
eight 50 amp-turn current loops of 1 m diameter. The 
loops are aligned to within 10 mrad of the detector axis 
using a sensitive field probe and an AC lock technique. 
Additional coils canceled the transverse components of 
the Earth's field and local gradients of 7.5 f/T/m. 

The vacuum chamber begins at the spin flipper with 
two meters of Be-coated flight tubes, through which the 
neutrons travel toward the collimator region. Two colli- 
mators of 6 cm and 5 cm diameter openings separated by 
2 m define the beam. These and 5 additional "scrapers" 
between them consist of rings of 6 LiF which absorb the 
neutrons. Behind each ring is a thick ring of high-purity 
lead which absorbs the gamma-rays from the reactor and 
those produced by neutron captures upstream. Between 
scrapers, the walls of the beam tube are lined with 6 Li- 
loaded glass to absorb stray neutrons. 
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FIG. 3. Two sheets of current-carrying wires create a mag- 
netic field of opposite orientation on each side. The field ori- 
entation changes so rapidly that the spin of a neutron passing 
through the current sheets cannot follow the field reversal, 
and the neutron polarization is reversed with respect to the 
magnetic guide field. Downstream the magnetic field and po- 
larization are rotated adiabatically from transverse in orien- 
tation to longitudinal. 



FIG. 4. Neutron beam intensity profile at the entrance to 
the detection region obtained from a scan across the beam 
face. 



B. Detector System 

Eight detectors surround the beam, each 10 cm from 
the beam axis as shown in Figure ^[ The octagonal ge- 
ometry places electron and proton detectors at relative 
angles of 45 and 135 degrees. Coincidences are counted 
between detectors at relative angles of 135 degrees. 



A fission chamber mounted behind a sheet of 6 Li-glass 
with a 1 mm pinhole aperture was scanned across the 
beam to obtain a cross-sectional profile of the intensity 
as shown in Figure ||. The neutron intensity was mea- 
sured before and after the experiment. To determine the 
polarization at the entrance to the detector, the beam 
passed through a second, single-sided, analyzing super- 
mirror directly in front of the scanning detector, and the 
ratio of intensities with the spin flipped and unflippcd 
was measured. The resulting fllppng ratio measures a 
combination of the neutron-spin-dependent transmission 
efficiencies of the two supermirrors and the neutron spin 
flipping eficiency. From this, and assumptions about the 
spin flipping efficiency, we can determine the product of 
polarization efficiencies for the two supermirrors (polar- 
izer and analyzer). When the upper limit of 100% spin 
flipping efficiency is used, a lower limit of the neutron 
beam polarization of 93% (95% CL) is found. This lower 
limit also includes the assumption that the flipping ratio 
for a pair of supermirrors identical to our analyzer would 
be less than that of a pair of supermirrors identical to 
our polarizer by a factor of 2 ± 0.5% pi] ]. 

Downstream of the detection region the vacuum cham- 
ber diameter increases to 40.6 cm, terminating with a 
6 Li-glass beam stop 2.8 m from the end of the detector. 
A 1 mm diameter pinhole at the center of the beamstop 
allows about 1% of the beam to pass through a silicon 
window into a fission chamber detector that continuously 
monitors the neutron flux. 



1. Electron Detectors 

The electron detectors are slabs (8.4 cm x 50 cm x 
0.64 cm) of BC408 plastic scintillator connected on each 
end to curved lucite light-guides that channel the light to 
Burle 8850 photomultiplier tubes. Each photomultiplicr 
tube is surrounded by a mu-metal magnetic shield and 
a pair of nested solenoids acting as an active magnetic 
shield. This combination of active and passive magnetic 
shielding had a factor of 10 less impact (0.5 /zT) on the 
guide field at the beam center than the mu-mctal alone. 

The scintillator thickness of 0.64 cm is just greater than 
that necessary to stop the most energetic (782 keV) of the 
electrons from neutron. The scintillators are wrapped 
with aluminized mylar and aluminum foil to prevent 
charging and to shield the detectors from x-rays and field- 
emission electrons in the vacuum chamber. For each seg- 
ment, the energy response was calibrated with cosmic-ray 
muons and conversion electrons from 207 Bi and 113 Sn (see 
Figure |(|) 

2. Proton Detectors 

Each proton detector has an array of 12 PIN diodes 
of 500 micrometer thickness arranged in two rows of 6. 
The diodes are held within a stainless steel high voltage 
electrode. Over each diode an open cylinder protrudes 
from the face of the electrode, shaping the field to focus 
and accelerate the protons as shown in Figure |. Thus 
each diode collects protons focused from a region 4 cm 
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FIG. 5. Basic detector geometry - an octagonal array of four each proton and electron detectors. 
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FIG. 6. Spectra produced during the energy calibration of 
the electron detectors. Shown are histograms of the charge 
collected in an individual phototube (PMT A) and of the total 
charge in the analog hardware sum of the two phototube sig- 
nal. The spectrum in PMT B (not shown) is roughly identical 
to that in PMT A. For the 113 Sn, a level-crossing discrimi- 
nator triggers on the analog sum signal. The peak visible is 
the 364 keV conversion electron. To suppress the contribu- 
tion of gamma-rays in the 207 Bi spectra the data acquisition 
is triggered on a thin scintillator placed between the source 
and the detector. The largest feature in these spectra falls at 
an energy of 882 keV (the 976 keV conversion electron energy 
minus the energy loss in the thin scintillator.) 
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FIG. 7. Geometry of the electrodes that accelerate and fo- 
cus the protons onto a PIN diode. 



x 4 cm even though it has an active area of only 1.8 cm 
x 1.8 cm. The diodes and their electronics are held at 
-30 to -40 kV. Between the electrode and the beam is a 
frame strung with 80 0.08-mm gold-plated tungsten wires 
that define a plane of electrical ground. Protons drift in 
a field-free region until they pass this plane, and then 
are accelerated by the high voltage and focused onto the 
nearest PIN diode below. Near both ends of the detector 
array are two cryopanels held at liquid nitrogen tempera- 
ture. Water vapor, released predominantly by the scintil- 
lators and other plastic components, is pumped onto the 
cryopanels to prevent condensation on the cooled PIN 
diodes. 

The charge in the PIN diode produced by each proton 
is amplified by 10 V/pC with a preamplifier mounted di- 
rectly behind the PIN diode. These circuits and the PIN 
diodes are cooled with liquid nitrogen to about 0°C to de- 
crease electronic noise. Preamplifier signals are processed 
in a custom VME-format shaper/ADC board with pro- 
grammable gain and operating mode parameters. The 
PIN diodes were calibrated with x-rays from an Am 
source as shown in Figure ^. 
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FIG. 8. Energy calibration spectrum of a PIN diode detec- 
tor using an 241 Am source. The FWHM is 2.9 keV at the 
59.5 keV line. 



3. Background 
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FIG. 9. Data acquisition components. 



The background in the detectors was primarily related 
to the beam or to the high voltage bias. Closing the beam 
shutter upstream of the neutron filter stops virtually all 
neutrons and about 1 /3 of the gamma-rays coming from 
the reactor along the beamline. With the shutter closed, 
the rates in each detector were less than 100 Hz, primar- 
ily from dark current, reactor gamma-rays, and cosmic 
rays. With the shutter open, the detectors see an in- 
creased gamma-ray flux primarily from neutron captures 
in the apparatus, triggering the detectors at less than 1 
kHz per electron detector and less than 1 kHz for all PIN 
diodes combined. This results in deadtimc less than 3% 
for the beam-related background. At its worst, the high- 
voltage-related background, consisting of x-rays, light, 
electrons, and ions, led to rates in the hundreds of kHz 
in the detectors. It was reduced at times by conditioning 
and cleaning of electrodes but varied by orders of mag- 
nitude during the run. 



C. Data Acquisition 

A block diagram of the data acquisition system is 
shown in Figure ^. The identification of neutron decay 
events is simplified by the fact that the proton signal 
is observed 0.5 /is to 2 /is after the electron signal. The 
recoil protons, with maximum energies of only 750 eV, re- 
quire this time to drift from the point of decay to the face 
of the proton detector. Events are accepted by the co- 
incidence trigger when the electron signal arrives within 
a coincidence time window ±r CO j nc /2 of a proton signal. 
The duration T CO i nc of this window was originally 14 /is 
and was shortened to 7 /is midway through the experi- 
ment to reduce the system deadtime. Each stored event 
contains the location (PIN diode) and energy for the 
proton event, location (electron detector), energy of the 
electron event, relative time between individual signals 



from the two phototubes in the electron detector, rela- 
tive time of arrival of the proton and electron signals, and 
the orientation of the neutron polarization. Every 30 sec- 
onds during the data collection, information is recorded 
from the system monitors which include system livetime, 
magnet currents, neutron flux at the beam stop, vacuum 
pressure, proton detector high voltage, and high voltage 
leakage current. Periodically, the data acquisition col- 
lects singles spectra from all of the individual detectors. 



III. EXPERIMENTAL RUN 

A. Data Collection 

The experiment was installed at the NCNR during De- 
cember 1996 and January 1997. From February through 
August 1997, 50 GB of data were collected and stored. 
The data are divided into 626 files representing contin- 
uous runs, typically four hours in duration. These are 
grouped into 125 series, within which running conditions 
varied little. For one week in August a systematic test 
was run in which the beam was distorted and the polar- 
ization guide field direction changed. The purpose and 



results of this test will be described in Section IV 



Instabilities in the proton detector high voltage made 
it impossible to operate all channels of the detectors at 
all times. Sometimes the electrodes simply would not 
hold the necessary voltage, and at other times a large 
spark or series of sparks would damage the electronics 
held at high voltage. Less than half of the data were 
collected when all four proton detectors were function- 
ing. Another limitation to the detector uniformity were 
variations in the measured proton energy deposited in 
the PIN diodes. In preliminary tests, the surface dead- 
layers of the PINs were measured to be 20±2 /ig/cm 2 as 
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FIG. 10. Energy spectrum in PIN diode III14, near which is 
mounted a weak 119 Sn source producing a 24 keV x-ray. The 
protons, accelerated to 36 keV but measured at less than 20 
keV, are visible between the background and the x-ray peak. 
The peak on the far right from a low-rate pulser input directly 
into the preamplifier is used to monitor gain and resolution. 
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specified by the manufacturer, Hamamatsu. In a dead- 
layer of this thickness a 35 keV proton loses 10 keV of 
energy. The proton energies measured during the experi- 
ment, however, were 12-18 keV, an average of 20 keV be- 
low the energy imparted to them by acceleration through 
34-38 kV (see Figure [u|) With widths (FWHM) of ap- 
proximately 10 keV, these peaks are not well separated 
from the background. High background rates necessi- 
tated the setting of thresholds at levels such that some 
neutron decay events were also rejected. This and the 
data acquisition deadtime were the primary limitations 
to the statistics of the experiment. A deadtime per event 
of 2 ms was necessary for stability of the system. Even 
with the reduction in length of the coincidence window, 
the high rate of background kept the system at 40-60% 
deadtime for most of the data collection period. 



B. Event Selection 

Figure [ll] shows an example of the relative time spec- 
trum for the coincidence data. The large center spike, 
originating mainly from multiple gamma rays produced 
by neutron captures in the apparatus, defines zero time 
difference. The neutron decay events are accepted within 
a window 0.35 /is to 0.9 /is after the prompt peak. This 
window contains the majority of the neutron decay pro- 
tons, while excluding the tail of the prompt peak and the 
low-signal-to-background tail of the proton peak. The 
background to be subtracted from these events is esti- 
mated using the rates in regions to either side of the decay 
and zero-time peaks. Events are also selected on the basis 
of measured proton energy to reduce the amount of back- 
ground to be subtracted. The energy range accepted is 
chosen solely by minimizing the fractional statistical un- 
certainty in the number of neutron decay events for each 
PIN diode-electron detector pair. Specifically, if -/Va is 
the number of coincidences counted by subtracting the 
background from the coincidences in the 0.35 to 0.9 /is 
window, the energy range is chosen to minimize 



gar* „ v 71 + V7 



(3.1) 



where / is the signal to background ratio in this energy 
range. This increases the overall signal to background on 
the 15 million good events from 0.8 to 2.5. 



IV. DATA ANALYSIS AND UNCERTAINTY 
ESTIMATION 

A. Determination of D from Coincidence Events 

For each PIN diode-electron detector pair in a given 
data series, the count rate can be expressed as 



N^ 1 = N e a e l [Kf + aK™ ± Pa ■ 
(AKf + BK% + DKfl)] . 



(4.1) 



where Nq is a constant proportional to the beam flux, the 
e a and e 1 are detector efficiencies for a PIN diode and 
electron detector respectively. The average of the neu- 
tron polarization vector over the detector volume, given 
by Prj,is assumed to be uniform and constant over time, 
lying along the direction of the 0.5 mT guide field. The 
± signs correspond to the two signs of the polarization. 
The factors K" 1 an d K" 1 are geometric factors derived 
from Equation 1.1 by integrating 1 and p e • p u /E e E u , 



respectively, over the beta-decay phase space, the neu- 
tron beam volume, and the acceptance of each electron- 
detector-PIN-diode detector pair. Similarly, the factors 
K™, Kg*, and K™, are obtained by integrating the vec- 
tors: p e /E e , p v /E v , and (p e x p v )/E e E v . 

We produce the following efficiency-independent asym- 
metries 



N1 



N1 



w 



From Equation |4.1| we get 



Pa ■ (AK 



N ai 



(4.2) 



(4.3) 
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PIN a PINb 

FIG. 12. The data from two PINs at the same z-position in 
a proton segment can be used to cancel the effects due to the 
electron and neutrino asymmetries. The coincidences shown 
by solid lines (EiPIN a and E2PIN(,) have approximately the 
same angle, a little less than 135°. These are referred to as 
"small-angle" coincidences. The "large-angle" coincidences 
for this pair of PINs (EiPINj, and E 2 PIN a ) are the dashed 
lines. 



where we use the definitions 



Kf 



etc. 



(4.4) 



Consider the two detector pairings PIN a -Ei and PINf,- 
E2 indicated in Figure [l2|. The corresponding values 
of K^ have opposite sign while and K™ have the 
same sign. We therefore combine asymmetries from two 
proton-electron detector pairings to produce the combi- 
nation 



,.62:al 



-Pa 



1 

" 2 



62 all 
W — W 



A(K^ - K^ 1 ) + B(K b i - Kg 



(4.5) 



(4.6) 



For uniform detection efficiency the difference (Kg — 
Kg) lies along the detector axis, z, while the differences 
(K^ 2 



detector axis 



— Kg) lie perpendicular to the 
For a polarized neutron beam with per- 



K^ 1 ) and (Kg 



feet cylindrical symmetry aligned with the detector axis, 



Kg 



= -K? 



and 



,.62:ol 



= PDK?? • z = -PDK°, 



(4.7) 



Departures from perfect symmetry and perfect alignment 
of the neutron polarization require that the A and B cor- 
relation terms be retained in Equation 4.6. The resulting 



systematic effects are discussed in Section [V C 

Additionally, as shown in Figure [l^, there are two 
classes of electron-PIN pairs: those that make an angle 
smaller than 135° (62 : al) or an angle larger than 135° 
(a2 : 61). We thus separate our data into a small-angle 
group and a large-angle group giving two statistically in- 
dependent results for each PIN-diode-electron-detector 
pairing. 



B. Monte Carlo Methods 

We use two Monte Carlo calculations to determine 
the values of K\, K a , Ka, Kb, and Kn. The results 
from these two completely independent calculations are 
in excellent agreement. In both calculations neutron de- 
cay events are generated randomly within a trapezoid- 
cylindrical geometry (i.e. a tube with divergence) that 
can be offset with respect to the detector axis. A realistic 
beam profile, representative of Figure ^, can be modeled 
by combining results from several different trapezoids. 
In one of the Monte Carlo calculations the tracking of 
protons and electrons is done with the CERN Library 
GEANT3 Monte Carlo package Q, while in the other 
tracking is implemented within the code itself. In both, 
the emiT detector geometry is specified with uniform ef- 
ficiency over the active area of each scintillator and over 
the square focusing region of each PI N d iode. 

The constants defined in Equation 4.1 are given by 



K° 



5 m X 



(4.8) 



where X = 1, p e -p v /E e E u , p e /E e , p u /E Ul and 
Pe x Pv/EeEv for x — 1, a, A, B, and D, respectively. 
We have studied systematic uncertainties associated with 
potential non-uniformities in the beta efficiencies and in- 
cluded them in the final uncertainty for the constants 
K" 1 . These constants (a total of 11, taking into account 
the three directions for each vector) are accumulated in a 
file that is read to calculate the factors v (Equation [0]) 
for different orientations of the polarization. 

Values of |K^-z| are used directly in the interpretation 
of the result for D. Variations among the PIN diode pairs 
of individual values of Kg 1 within a given proton segment 
are negligible, and average values (\Kd ■ z\) can be used. 
They are found to be 0.424 ± 0.010 and 0.335 ± 0.020, 
for the small- and large-angle coincidences, respectively. 
The uncertainties are primarily from uncertainties in the 
geometry of the beam. Values for the other K£ l are used 
in the estimation of systematic uncertainties described in 
the following section. 



C. Discussion of Systematic Uncertainties 

The largest of the systematic effects c an b e shown to 
be the contributions to the v (Equation [4.9| ) that arise 
due to the misalignment of the neutron polarization with 
respect to the detector axis. A transverse component 
of the polarization produces a significant contribution 



to v 



b2:al 



because the vector differences K 



1,2 



K^ 1 and 



K^ — K^ 1 are predominantly perpendicular to the de- 
tector axis. (For example, K^ — K^ 1 is proportional to 
the integral of p e (E\) — p e (E2) and is directed horizon- 
tally to the left in Figure pi The difference Kg - K^ 1 is 



antiparallel to K 



^ A 



For an azimuthally symmetric 
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neutron beam, it can be shown that for each proton de- 
tector segment (labeled with subscripts 77= I, II, III, IV) 
the weighted average of the v av -Pi for all large or small 
detector pairs can be expressed as 



J/s _ 



= PD(K l l, s ■ a) + a\' s sin 6^ sm(<j> v - 



l/s, 



(4.9) 



where 9 a and <p a are the polar and azimuthal angles of rj, 
and r ,=O°, 90°, 180°, and 270° respectively for detectors 
I, II, III, and IV. This dependence can be derived ana- 
lytically for zero beam radius and is confirmed by Monte 
Carlo simulations for symmetric beams of finite radius. 
The coefficients a n measure the combined effects of the 
A and B correlations for each proton detector segment. 

If the symmetry of the four sets of proton detectors 
were perfect, i.e. a\ — a\\ = «in = arv, the contribu- 
tions due to the A and B coefficients would average to 
zero, and Equation 4.7 would be valid, even with a po- 
larization misalignment. In the absence of perfect sym- 
metry, these contributions do not cancel when the four 
proton detectors are combined, and a false D contribu- 



tion would result from the application of Equation 4.7 



This false D is proportional to the product of two effects 
that are both small: the misalignment of the neutron po- 
larization with respect to the detector axis (9 a ) and the 
departure from perfect symmetry of the proton detectors 
(Aa = l/2(ai — am) + l/2(an — arv))- Such an effect 
is called the "tilting asymmetric transverse polarization" 
effect, or "Tilt ATP" 

The ATP effect was intentionally amplified for a sys- 
tematic test, run with transverse polarization (6 a = 90° 
4>a = 4>iv = 270°) and a distorted neutron beam. The 
neutron beam was distorted by blocking half of the beam 
with a neutron absorber placed upstream near the spin 
flipper. The results of this test are shown in Figure 
This demonstration that the experiment can measure an 
asymmetry consistent with the Monte Carlo calculation 
serves as a strong check on both the operation of the 
detector and the validity of the analysis method. 

A false D also arises if the polarization has transverse 
components not described by a simple tilt. The form of 
Equation 4.9 shows that a net azimuthal component of 
a also results in a contribution to that does not av- 
erage to zero when data from proton segments I-IV are 
combined. This effect, referred to as a "twisting asym- 
metric transverse polarization" ( "Twist ATP" ) is shown 
by Monte Carlo simulations to be less than 10~ 4 for az- 
imuthal polarizations of less than 1 mrad. For this rea- 
son, all sources of guide field distortion are kept to less 
than 1 mrad, and materials of low magnetic permeability 
(less than 0.005 /io) were used in the detection region. 
There are exceptions to this requirement, however the 
net effect of all additional permeability was measured to 
produce less than 1 mrad of distortion of the guide field 
anywhere in the detector region. 

Variations in the neutron flux ($) and polarization (P) 
that depend on neutron helicity yield a false D. For this 
experiment the effects due to misalignment of the neutron 




-1.0 



small angle coincidences: 

• systematic test data 

O simulation 
large angle coincidences: 

■ systematic test data 

□ simulation 



T 



1 1 1 1 

50 100 150 200 250 

azimuthal angle <t> of proton detector (degrees) 



FIG. 13. During the systematic test with a 90° polariza- 
tion tilt, the false D in each proton segment is clearly vis- 
ible. (The proton detector segment IV at 0rv = 270° was 
not operational.) Also shown are the results of Monte Carlo 
simulations of the systematic test conditions. The error bars 
shown are purely statistical, and are not an accurate estimate 
of the total uncertainty in the calculation, of which the largest 
contribution is uncertainty about the shape of the beam. 



spin are small so that these systematic effects, to 
first order in A<i>/<I> and AP are 



D /aise (A$) = -^PDo ■ {A(K A ) + B(K B )), (4.10) 



and 



D false (AP) = APD& ■ (A(K A ) + B(K B )). (4.11) 

Here A$ = $ T -$j , and AP = P + -P_ . (P in Equations 
16] and [O] would be replaced by P = (1/2) [P+ + P_].) 
The (Ka) and (Kb) are average values for all PIN-diode- 
electron-detector pairings. 

Our data provide an upper limit of 0.002 for Pa ■ 
(A(Ka) + B(Kb))- We combine this with neutron 
flux monitor data for < 0.004, concluding that 

Df a i se (A$) < 8 x 10 _6 P. The flipping ratio mea- 
surement has been used to derive a lower limit on the 
spin flipper efficiency of 82% so that AP < 0.2, and 
D false (AP) < 4 x 10~ 4 D. We conclude that both effects 
are negligible in this measurement. 



D. Results 



A final value of D 



l/s 



v v /(PKf 



z) is found sep- 



arately for large angle and small angle pairings of each 
proton segment. The quantities v v are the weighted av- 
erages of all PIN-electron detector pairs, v l ^ s (ai : 
within each proton detector segment. Use of the weighted 
averages is justified because the systematic uncertainties 
described above have negligible variations among the PIN 
diode pairs in a given detector. The individual proton 
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-5 ■ 



-15x1 CP 



• small angle 
coincidences 

■ large angle 
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"I 1 1 1 T 

50 100 150 200 250 

azimuthal angle $ of proton detector (degrees) 

FIG. 14. Results from each of the four proton segments 
for small-angle and large-angle PIN-electron detector pairs. 

■£) 



Dl/ S = Vr,/{PK% 



l/s 



Error bars are statistical. 



segment data (v^) are then combined in an arithmetic 
average so that the sinusoidal variation given in Equa- 
tion 



4.7 cancels to first order in misalignments, i.e. 



IV 



(4.12) 



The error for D l / S includes the uncertainty in the values 



of K% s ■ z. 



1 



IV 



' D l/s 



r/=I 



IK', 



IK 



(4.13) 



Data for each proton segment are displayed in Figure [IJ, 
where we plot values for the eight individual d\ { s . (The 
sinusoidal variation is predicted by Equation |4.9| and also 
seen in the test data of Figure [l3|, where the amplitude 
is 100 times larger.) 

The two independent measurements for small angle 
and large angle PIN-electron detector pairs can be com- 
bined in a weighted average. 



D 



D*/a 2 D . 



D l /a 2 



D< 



1/a 2 



(4.14) 



D l 



The full uncertainty includes the uncertainty from the 
average neutron beam polarization. 



1 



P 



(4.15) 



The data are also analyzed by breaking each series up 
into individual runs and combining PIN-electron detec- 
tor pairings in the same way. The results of these analy- 
ses are consistent. The final result is (— 0.6± 1.2) x 10 -3 , 
where we have assumed the neutron polarization is P = 



(96 ± 2)%. This is derived fro m our measurement of flip- 
ping ratio described in Section II A, with the assumption 
that the allowed range (93% < P < 100%) spans 2a p. 

Finally, we use the scaled results from the systematic 
test data (Figure |l3|) combined with Monte Carlo sim- 
ulation studies to estimate the uncertainty of the Tilt- 
ATP systematic effect. For the test data, proton detec- 
tor IV (0iv = 270°) was not operational. In calculating 
D for the test data, only values fro m de tectors I and 
III can therefore be used in Equation 4.12 with a result 
of i(L> T + D m ) = (-6.5 ± 1.4) x 10- 



Monte Carlo 

simulations show that for a beam of radius 3 cm, the 
sin(</) r( — ) behavior of Equation |4.9| is modified so that 
Dtest = \{D l + Dm)/ 1.6 = (-4.1 ± 0.9) x 



10~ 2 . This 

can be scaled by sint9 CT , the ratio of polarization mis- 
alignments for the data and test runs. The individual 
values of D l ^ s shown in Figure [li] are used to determine 
6 a = (9±3) x 10~ 3 radians for the data run. This provides 
an upper limit for the uncertainty on the Tilt- ATP sys- 
tematic effect of £>(Tilt ATP)< D test sint9 CT < 5.2 x 10~ 4 . 
Though we use the test results to estimate this false D 
effect, we expect the cancellation due to beam symme- 
try to be more complete for the data run because the 
test beam was intentionally distorted. We therefore con- 
sider this upper limit to be a conservative estimate of the 
largest possible false D effect [^6| The contributions to 
the statistical and systematic uncertainties are given in 
Table g 



Sources of uncertainty 


Contribution (xlO 4 ) 


Statistics 


12 




Tilt ATP 


5 




Twist ATP 


<1 




Flux variations 


neglif 


;ible 


Polarization variations 


neglif 


;ible 



TABLE II. Contributions to the uncertainty. 



V. SUMMARY AND CONCLUSIONS 

The apparatus used to perform a measurement of the 
D-coefficient in the beta-decay of polarized neutrons has 
been described. The data using the emiT detector have 
been analyzed using a technique that is insensitive to 
the nonuniform detection efficiency over the proton de- 
tectors. The initial run produced a statistically limited 
result of D = [-0.6 ± 1.2(stat) ± 0.5(syst)] x 10" 3 . This 
result can be combined with earlier measurements to pro- 
duce a new world average for the neutron D coefficient of 
—5.5 ± 9.5 x 10 -4 , which constrains the phase of gA/gv 
to 180.073° ± 0.12°. This represents a 33% improvement 
(95% C.L.) over limits set by the current world average, 
and correspondingly further constrains standard model 
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extensions with leptoquarks p8j. The result is also in- 
teresting in light of upper limits provided by the neutron 
and 199 Hg electric dipole moments on T-odd, P-even in- 
teractions such as left-right symmetric models and exotic 
fcrmion models. 

A second run is being planned with strategies to im- 
prove the statistical limitations related to background ex- 
perienced in the first run. Our study of systematic effects 
presented here shows that the largest is the tilt- ATP ef- 
fect. The uncertainty on this effect can be reduced signif- 
icantly with more data taken i n the transverse polariza- 
tion mode described in Section IV C| . With the planned 
improvements in place, it will be feasible to improve the 
sensitivity to D to 3 x 1CK 4 or 
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